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ABSTRACT 

We made use of the two latest sets of rotational measures (RMs) of extra-galactic radio sources, 
namely the NRAO VLA Sky Survey rotation measures catalog, and a compilation by Kronberg and 
Newton-McGee, to infer the global structure of the Galactic magnetic field (GMF). We have checked 
that these two data sets are mutually consistent. Given the existence of clear patterns in all-sky RM 
distribution we considered GMF models consisting of two components: disk (spiral or ring) and halo. 
The parameters of these components were determined by fitting different model field geometries to 
the observed RMs. We found that the model consisting of a symmetric (with respect to the Galactic 
plane) spiral disk and anti-symmetric halo fits the data best and reproduces the observed distribution 
of RMs over the sky very well. We confirm that ring disk models are disfavored. Our results favor 
small pitch angles around ^ — 5° and an increased vertical scale of electron distribution, in agreement 
with some recent studies. Based on our fits, we select two benchmark models suitable for studies of 
cosmic ray propagation, including the ultra-high energies. 

Subject headings: Galaxy: structure -ISM: magnetic fields- methods: data analysis 



1. INTRODUCTION 

A realistic model of the Galactic magnetic field (GMF) 
is needed for various applications, such as direct dark 
matter searches, studies of cosmic rays, and others. In 
particular, the GMF plays a crucial role in the propaga- 
tion of ultra-high energy cosmic rays (UHECRs) — those 
with energies in excess of 10^^ eV. Their deflections (as- 
suming these are charged particles) in the GMF are large 
enough to prevent the identification of some, perhaps 
most sources directly from the UHECR data. However, 
in the case of protons these deflections are sufficiently 
small to be corrected for, provided the GMF is known 
with enough accuracy. Thus, the possibility of astron- 
omy with charged particles at ultra-high energies may 
depend crucially on our knowledge of the GMF. 

The magnetic field (MF) of our Galaxy is thought to 
conta i n both regu lar and turbulent components (see, e.g., 
iBeckl ()2001l . l2008f) ). While the regular component is sub- 
dominant in strength by a factor of a few, it is expected to 
give a dominant contribution in integral quantities such 
as total deflections of ultra-high energy particles. Unlike 
the turbulent component, the regul ar component adds up 
coher ently and comes to dominate (|Tinvakov fc TkachevI 
l2005f ) when integrated over many coherence lengths. Es- 
timating the strength and global structure of the regular 
MF of our Galaxy is the main purpose of this paper. Our 
primary interest is on MF outside the thin Galactic disk, 
since that is the field that determines the deflections for 
most arrival directions of UHECRs. 

The regular component of the GMF ha s been studied 
previously e.g. (fS iniard-Normandin & Kronberg 1983 
ISofue fc FuiimotorflQSa IHan fc Oiad 119941: iFrick'et ahl 



200lt I Vailed [2005t IHan et all IpOOl: I Brown et all [20071: 
Sun et al .1 120081 : IJansson et al.l 1200(1 . and in particular 



in the context of UHECR propagation (jStanev' Il997t 
iTinvakov fc TkachevI [200l iProuza fc Smi da 2003). The 
advance of this paper is that it is based on improved 
data on rotation measures (RMs) and the ionized gas 
distribution. We also develop improved numerical model 
commensurate with this improved data. 

Observations of other spiral galaxies reveal MFs of a 
few ^lG strength, coherent over distances of severa l kilo- 
parse cs (see, e.g., (Beck 2008; Krausc 2009; Brau n et al.l 
120101 )). Observations from the face-on and edge-on per- 
spectives indicate that these fields consist of at least two 
components: the 'disk' one which is concentrated to- 
wards the galactic plane and has largely a spiral pattern, 
and the 'halo' one that resides at some distance from the 
disk plane and has radial size comparable with the size 
of the galaxy. As the Milky Way seems to be a typical 
member of the entire class of spiral galaxies, MFs of sim- 
ilar strength and configuration should be expected in our 
Galaxy. 

There are several methods to estimate the MF 
in the Ga laxy, suc h as Zeeman splitting mea- 
surements (iCrutched Il999f). i nfrared polarization 
studies (jNishivama et al.l [20100 . synchrotron po- 
larization and intensity surveys fWiclcbinski 200^ 
'Sun et al.' '2008^; 'Janss on_et all [2009; Ja ffe et al.i 12010: 
Ruiz-Granados et al. 2O1O0, starlight polarization stud- 
ies (Heiles 1996) and, finally, observation of the Faraday 
RMs of different Galactic Jpulsars) an d extragalactic 



radio sour c es (e. g . r^s 

Han et al.' 120061: [Men et al 



(iSim ard-Norman din et aLl 
_ ■ _ ■ ■ , et al.l i20q8i: iMao et al.l 

Kronberg fc Newton-McGed l2011f) . Except for 
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Zeeman measurements all these methods are somewhat 
indirect as they involve additional assumptions such as 
density of warm ionized medium or concentrations of 
cosmic ray electrons. In this paper we use recent data 
on the Faraday RMs. 

When propagating through a magnetized plasma, the 
polarization plane of a linearly polarized electromagnetic 
wave of wavelength A rotates by the angle Ai/j propor- 
tional to the square of the wavelength, 



(1) 



Multi-frequency observations thus allow one to measure 
the value of RM. The RM can be expressed in terms 
of the properties of the intergalactic/interstellar medium 
(IGM/ISM) and permeating MFs: 



RM = 0.81 



npBiidl 



(2) 



where n, 
cni^-^, B 



is the density of free electrons measured in 
is the component of the MF parallel to the line 
of sight measured in fiG (positive when directed towards 
the observer), and D is the distance from the observer 
to the source in pc. Combined with data on the Galactic 
free electron distribution. Equations ([T]) and ([2|) can be 
used to estimate the GMF from the Faraday RMs of a 
large number of sources. 

Two types of sources are normally used to measure 
Galactic Faraday rotation: Galactic pulsars and extra- 
galactic radio sources. Pulsars have an advantage that 
their intrinsic RM values are negligibly small. However, 
the use of pulsars requires the knowledge of their dis- 
tances (usually inferred from their dispersion measures 
DM), which introduces additional assumptions about 
the ISM density distribution. In addition, pulsars are 
concentrated around the Galactic plane, and their total 
number (~ 500) is currently much smaller than the num- 
ber of extragalactic radio sources with measured RMs. 
For these reasons in this paper, we concentrate on the 
RMs of the extragalactic radio sources. 

Our general strategy is the following. We consider sev- 
eral types of analytical GMF models, each containing a 
finite number of parameters. For each model we calculate 
the expected RMs, and fit them to the observations by 
the binned method. The quality of the fit is then used 
to determine the best values of the GMF parameters and 
their allowed ranges. 

The paper is organized as follows. Section [2] describes 
the data. Sections [3] and S] describe the implemented 
models and our method. The results are presented in 
Section[51 Section [5] contains the results and conclusions. 

2. DATA 

In the present work we use two sets of the RM data. 
T he first one is the com pilation of RMs recently obtained 
bv lTavlor et"all H2OO90 who reanalyzed NRAO VLA Sky 
Survey (NVSS) data. NVSS is the largest by number 
survey of p olarized radio sources at declinations > —40° 
(jCondon et al. 1998). It was made in two nearby bands, 
1364.9 and 1435.1 MHz; each having a width of 42 MHz. 
Simultaneous observations in these two different bands 
give estimates of the RMs of the sources. The total num- 
ber of the observed sources was 37,543 with the mean 
error, estimated by the authors, of ~ 11 rad m~^. The 
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Fig. 1. — Locations and RM signs of the NVSS sources. Red 
(blue) color corresponds to positive (negative) values of RM. 
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Fig. 2. — Upper panel: locations and RM signs of 2247 sources 
from KNMll. Red (blue) color corresponds to positive (negative) 
values of RM. Lower panel; the same data smoothed with a re so- 
lution of 21° (adopted from lKronberg Newton-McGei l|20TlT) 1. 

sky map of the NVS^l set is shown in Figure [1] which 
displays the location and RM sign of each source. The 
data have a large blind spot corresponding to declina- 
tions below ^ -40°. 

An advantage of the NVSS RM catalog is the large 
number of sources. However, observations being made 
in two close bands lead to relatively large individual RM 
errors. Also because in general RM is not always con- 
stant in A^, it can have modest variations over a wider 
range of A. Thus an NVSS RM, measured in just one 
small wavelength range at 21 cm (A = 20.9 - 22.0 cm), 
can differ from an RM measured over a wider wavelength 

^ We henceforth refer to NVSS rotation measure catalog as just 
NVSS. 
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range which tends to average out fluctuations in RM(A^). 
Such systematic differences will add to the mean error of 
~ 11 rad m~^ quoted above by up to ~ 13 rad m~^, and 
will contribute to the width s of the RM comparison plots 
with the results from the IKronberg fc Newton-McGed 
(|2011[) compilation (KNMll) shown in Figure |3] below. 
On the other hand, because A2 — is small, ambiguities 
of n TT in the NVSS RM's are unlikely to occur, especially 
at |6| > 10°. 

It can be easily seen that the NVSS data set should be 
treated with care when considering the sources with large 
observed RMs, which mainly reside near the Galactic 
plane and in the region of inner Galaxy. In order to avoid 
this problem, we use the NVSS data only at Galactic 
latitudes \b\ > 10°. 

The second all-sky d ata set we used is the compilation 
by IKronberg fc Newton- McGed ()201lD (2257 sources). It 
consists of ~ 1500 revised and statistically more accurate 
extragalactic radio source RMs of P. P. Kronberg et al. 
(in preparation). These are derived from polarization 
measurements made over a large "baseline" of A^. The 
remaining RM 's are taken from the Canadian Galactic 
Plane Survey dBrown et al.ll2003D. th e Southern Galac- 
tic Plane Survey dBrown et al.l 20071) . and two smaller 
published li sts from iKlein et all (|200l (108 RMs) and 
i)Mao et al.l I2008D (68 RMs). The extragalactic radio 
sources in the KNMll set were carefully observed at mul- 
tiple frequencies and give RMs largely free from ambigu- 
ities. Errors of the RM measurements in the KNMll set 
are typically of order ~ 3 radm~^. We use this set for 
control. 

The sky map of the KNMll RMs is shown in Fig- 
ure[2l upper panel. The lower panel shows the same data 
smoothed over the circle of 2 1° centered at the location of 
each s ource as described by IKronberg fc Newton-McGe3 
(Unni). One may observe a generally good agreement in 
the RM sky distributions of the NVSS and KNMll. This 
is particularly visible when comparing Figure [T] with the 
lower panel of Figure [2] 

The RM sky structure is rather peculiar. The re- 
gions of positive and negative RMs are roughly sym- 
metric with respect to the Galactic plane in the outer 
Galaxy 90° < I < 270° and anti-symmetric in the in- 
ner Galaxy. Note, however, that the boundary between 
positive and negative RMs in the inner Galaxy is dis- 
placed away from the Galactic plane (toward negative b 
for < Z < 90°), that is, RMs do not change sign across 
the Galactic plane itself. Th i s fact was observed, e.g., in 
IKronberg &: Newton-McGed ()20lH) . This overall struc- 
ture of RMs suggests, qualitatively, a Galactic field struc- 
ture that consists of two components, with the disk field 
being symmetric in b with respect to the Galactic plane 
and the halo field anti-symmetric. 

Although the second data set is significantly smaller, 
it has much smaller individual RM errors and can be 
used in two ways. First, it allows for an independent 
evaluation of the quality of the NVSS data. In order 
to assess the accuracy of the NVSS RMs we identified, 
by positional cross-correlation, the sources that are com- 
mon to both catalogs. We did it separately for the two 
regions |&| > 10° and |&| < 10°. We found 1338 pairs 
and 306 pairs of sources, respectively, separated by less 
than 0°.05. For these sources, we built distributions of 




1 

f 




1 






lbl<10° 


J 







-1000 



-500 





A RM 



Fig. 3. — Distribution of differences of RMs as given by the NVSS 
catalog and by KNMll. The upper curve, for |6| > 10°, shows a 
Gaussian fit for the range (-200, 200) rad ra~^. The lower curve, 
for |fe| < 10°, shows the corresponding result for the full range 
(-1000, 1000) rad m'^. 

differences between RMs given in the NVSS catalog and 
in the set by KNMll. The difference distributions are 
shown in Figure [3] 

In both regions, the distributions of ARM have peaks 
centered at zero and well fitted by Gaussians with dis- 
persions cr ~ 13 rad m^^ and ct ~ 36 rad m^^, respec- 
tively. In the region |6| > 10° the tails contain ^ 25% 
of points and are confined within ±100 radm^^. The 
effective error introduced by these tails is of the order of 
the Gaussian part. Thus, in this region the NVSS data 
can be used for our purposes. 

On the contrary, in the region |6| < 10° the tails (con- 
taining similar fraction of points) are spread over a very 
wide range ±900 radm~^. Such tails would introduce 
errors in the low-6 bins which are too large to constrain 
the fit. For this reason we exclude the bins with |6| < 10° 
from the NVSS fit, as will be described in Section 2] 

The KNMll RM set provides statistically more accu- 
rate RMs in the Galactic plane, including areas where 
the NVSS data do not give sky coverage. This allows us 
to check whether the models which fit the NVSS data 
also fit the RM data in the Galactic plane. 

3. GMF AND ELECTRON DENSITY MODELS 

We adopt a general model of the GMF consisting 
of two different components: a disk field and a halo 
field (|P rouza fc Smidal[200a ISun et al.l 120081) . Each of 
these two components is parameterized independently. 
According to the qualitative features of the RM distri- 
bution discussed above, for the most part of this work 
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Fig. 4. — Sketch of the structure of the galactic magnetic field. 
The disk field in this sketch has reversals of the MF in adjacent 
arms. It is shown that the halo field is antisymmetric with respect 
to the galactic plane: in the solar neighborhood the direction of 
the halo field is anti-coincident with the direction of the disk field 
above the galactic plane and coincident below. We assumed that 
the disk field is symmetric with respect to the galactic plane, in 
the sense that the sign of the disk field is locally preserved across 
6 = 0. 

we take the disk field to be symmetric with respect to 
the Galactic plane and the halo field anti-symmetric 
(other possibilities will be discussed in Section [S]). The 
general layout of the model is shown in Figure IH 



3.1. The disk field 

We consider several disk GMF models. The first 
one is the widely used logarithmic spiral model, 
e.g.. (iSimard-Normandin fc Kronber^llQSOHHan fc Qiao! 
[I99l IStanevI 119971 iTinvakov fc Tkachevri2002D . There 
are two versions of this model depending on whether the 
direction of the field in two different arms is the sam43 
(axisymmetric, or ASS model) or opposite (bisymmetric, 
or BSS model). They are defined as follows: 



Bg — Bcosp, Br — Bs'mp, 



B{r,e,z) = B{r) 



cos 6* — 6 In ■ 



Rr. 



exp(-|z|/zo), 
(3) 



or 



B{r,0,z) = B{r) cos 



bin — + exp(-|z|/zo), 



(4) 



61n 1 



for ASS and BSS models, respectively. Here 

d \ TT 

6=1/ tanp, 

where p is the pitch angle and d is the distance to the 
first field reversal. Negative d means that the nearest 
reversal occurs in the direction to the Galactic center, 
positive corresponds to the opposite direction. Rq is the 
distance to the Galactic Center (we adopt Rq = 8.5 kpc 
in this paper). The amplitude of the GMF is a function 
of the radial coordinate r. 



B{r) = 



Br 



Br 



Rr. 



r cos (j) ' 

Rq 



r>Rc 
r < Rr 



(5) 



Rc cos ( 

^ In the work of IJansson et aLl 120091 ) this model is dubbed as 
DSS- Disymmetric Spiral. 



where Rc is the radius of the central region where the 
disk field is assumed to have constant magnitude. In 
the course of the simulations, the parameters were varied 
within the following ranges: 



parameter 


min 


max 


Rc 


3 kpc 


6 kpc 




0.5 kpc 


1.5 kpc 


d 


-1.4 kpc 


1.4 kpc 


P 


-15° 


15° 



The second model of the disk field was the axisymmet- 
ric fi eld with reversal s in concentric rings. It was taken 
from lSun et all (l200l : 



B{r,9,z) = Di(r,z)D2(r) 
Rq H 



D2{r) = 




r>Rc 
r<Rc 



r > 7.5 kpc, 
6 kpc < r < 7.5 kpc, 
5 kpc < r < 6 kpc, 
r < 5 kpc. 



(6) 



(7) 



where Rq is radial scale of the MF. We studied the fol- 
lowing parameter space: 



parameter 


min 


max 


Rc 


4.5 kpc 


5.5 kpc 


Ro 


8 kpc 


15 kpc 


zo 


0.5 kpc 


1.5 kpc 



We also varied slightly the radial positions of the inner 
MF reversals (the boundaries of regions of different sign 

in Equation Q) in the ranges 4.8 5.2 kpc, 5.8 — 

—6.2 kpc and 7.3 8.0 kpc. 

The magnitude of the disk MF at the vicinity of the 
solar System Bq was taken to be 2 fiG. Note that, in 
the total GMF, the disk and halo components will be 
combined with arbitrary coefficients (to be determined 
by fitting), so the choice Bq = 2 /xG at this point is 
merely a matter of definition. 

3.2. The halo model 

The basic halo mod el was taken from (jProuza fc Smfdal 
[20M[S^menil[M)a ) 

B^ir,z) = 



B, 



H 



1 + 



R^ 



77y exp 1 - — g- , (8) 



R^ 



where the di rection of the field is reversed in the Souther n 
hemisphere ()Brandenburg et al.l 119921 : iHan et al.lll997D. 

We repeat here the following definitions from'Sun et alj 
(2008): B^ is the halo MF strength, i?^ is its radial 
scale, parameter Zq defines the vertical position of the 
halo, z^ is its vertical scale; this scale could differ in 
directions to and away from the galactic plane, \z\ < 



and \z\ > 



-,H 



respectively. We will refer to these 



parameters as z^-^-j and zf(2)' correspondingly. Here and 
below the index " H" denotes parameters belonging to the 
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ha lo. The original v alues of the parameters in the model 
of lSun et alJ (120081) were: ^ 10 /iG, R" ^ ^ kpc, 



-0 
■.H 



1.5 kpc and zf^ = 0.2 kpc for 



< Zq , otherwise 

^1 - 0.4 kpc ( z^^^ = 0.2 kpc and z'^^^^^ = 0.4 kpc). 
We also tested halo models with slightly different radial 
profiles: 

B^{r,z) = 



and 



H 



-,H 



exp 



B^ir,z) = 



B, 



exp 



r-Rff 



(9) 



Rt 



(10) 



We investigated the following parameter space: 



parameter 


min 


max 


Rq 

Rx 


3.5 kpc 
1 kpc 


15 kpc 
5 kpc 




0.8 kpc 
0.2 kpc 


3.5 kpc 
0.4 kpc 


^1(2) 


0.3 kpc 


0.5 kpc 



Given that the halo field has not been firmly established 
as entirely azimuthal, we further studied the basic halo 
model, varying the pitch angle within the range —10° < 



< 10° 



3.3. 



The electron density 

Calculation of the RMs, Equation ([2]), requires a 
knowledge of the Galactic electron density rie- In this 
paper we adop t ed th e NE2001 model developed by 
ICordes k, Laziol (j2002D . However, there is mounting 
evidence that the original vertical scale of this model 
is insufficient and pulsar observations are better de- 
scribed by the electro n density with larger vertical scale 
(|Gaensler et al.|[2008l) . We tried both the original [zh = 
0.95 kpc, Ue = 0.03 cm~^ at z = 0) and modified versions 
{zh = 1.8 kpc. He = 0.014 cm-3 at z = 0) of NE2001. 

The preliminary simulations with both versions of 
NE2001 were made with a coarse grid of GMF parame- 
ters. The original model performed systematically worse 
than the modified one, so we proceeded to the full sim- 
ulation with the modified version of NE2001. The re- 
sults presented below were obtained with this modified 
version. Apart from the test of two vertical scales as 
described in Section [5l we did not vary any other param- 
eters of NE2001 during the fits. 

4. METHOD 
4.1. Binning 

RMs of individual sources are affected by uncertain- 
ties due to the intrinsic source contributions to the RM 
and by small-scale fluctuations of the GMF. The effect 
of the coherent component of the GMF can be revealed 
by averaging over a large number of nearby sources. The 
angular bin size should be small enough so that vari- 
ations of the coherent component of GMF within the 




Fig. 5. — Binned celestial sphere. The bins marked in gray are 
excluded from the fit (see Appendix |A] for details). 

bin can be neglected, and it should include a sufficient 
number of sources to suppress source-to-source fluctua- 
tions. Applying these requirements, we chose the fol- 
lowing binning: the celestial sphere is divided into 18 
bands of 10° width, from -90° to +90° in Galactic lat- 
itude b. Each band is further divided into bins so that 
individual bins would have roughly equal area, about 100 
deg^, as shown in Figure [S] So, there are 36 bins in the 
9th band (0° < & < 10°), but only four in the flrst one 
(80° < b < 90°). The total number of bins is equal to 
422. After averaging, all the NVSS data are reduced to 
422 values of average RM in each bin. No outliers were 
removed at this stage. Figures [51 and [7] shows the av- 
eraged RMs in the eight bands with 10° < b < 50° and 
—50° < b < —10°, respectively. The coherent structure 
is clearly visible. 

The GMF models were treated in a similar way. For 
each model we simulated 20,000 valucfO of RMs at ran- 
domly chosen locations and then applied the same bin- 
ning and averaging procedure to these sets. The resulting 
averaged RMs were compared to observations. 

A somewhat different procedure was used for the set 
of sources from KNMll. To check the disk flcld, we used 
the sources residing within 1° from the galactic plane. As 
their number is more limited (about 240) and they are 
located irregularly within their respective bins (10 deg^ 
in this case), we simulated model RM values in their 
exact locations to avoid an additional error that would 
be produced by different distributions of simulated and 
real sources within the bin. 

4.2. Errors 

Important for this analysis is the weighting of different 
bins or, equivalently, the bin-by-bin estimation of errors. 
There are several sources of scatter of individual RMs 
within the bin: internal RM of the source, an anomalous 
extragalactic RM component, random component of the 
Galactic field, variation of the regular GMF over the bin, 
fluctuations of the Galactic electron density Ue, and the 
measurement error. All these factors contribute to de- 
viations of the observed average value of RM from that 
predicted by the GMF model, and determine the error 
that has to be assigned to the bin. 

If all the above factors were random, it would be pos- 
sible to estimate the error by statistical means. Note, 
however, that some of these factors (e.g., random com- 

^ We checked that increase in number of simulated sources from 
20,000 to 40,000 does not lead to any significant changes. 
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Fig. 6. — Averaged RMs in four bands 10° < 6 < 50°. The error 
bars show the errors used in the fit. 

poncnt of GMF with coherence length comparable to the 
bin size, fluctuations of n,,) may act coherently over the 
bin. These factors are much more difficult to treat rig- 
orously. We therefore take a phenomenological approach 
and assume that these factors produce an error in the 
bin that is some fraction of the root mean square (rms) 
of the RMs in that bin. In Appendix |X] we develop a toy 
model where such proportionality takes place. 
Taking account of this complication, we have estimated 



Fig. 7. — Averaged RMs in four bands 
error bars show the errors used in the fit. 



-50° < 6 < -10°. The 



the errors of RMs in the individual bins as follows. By 
considering the sources in the directions of north and 
south Galactic poles where the effect of the regular MF is 
expected to be suppressed, the dispersion in the observed 
RMs of quasars is estimated as Cqso — 15 rad-m~^ (see 
Appendix \K\. If the rms ai of RMs in the ith bin did 
not exceed this value, we set the error of average RMs to 



Cqso/v^Vi, where Ni is the number of sources in the 
bin. The reason is that in this case there is no indication 
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of coherence effects, so one may assume that the error 
is due to stochastic factors such as the internal scatter, 
measurement error, and short-scale random MFs. 

In the opposite case we assume that the variance in 
the observed RMs is (partly) due to the presence of some 
unknown features that could act coherently over the bin, 
but cannot be accurately modeled. We assigned the error 
of (Ti/3 to these bins, where the factor 1/3 is taken from 
the toy model of Appendix [X] We have checked that the 
best-fit parameters do not change significantly when this 
factor is varied in the range 1/2 — 1/4. This uncertainty 
in the absolute normalization of errors we cannot reduce 
further. 

After a bin-by-bin examination, we have increased the 
errors in bins where there are clearly visible anomalies in 
the observed data (see the list of these bins below in the 
Appendix |A|). These problems are: (1) existence of small 
cluster of sources with RMs sharply different from the 
surrounding values (2) a strongly irregular distribution 
of the observed sources across the bin (due to, e.g., deep 
surveys within small regions) that could bias the average 
compared to the case of a random distribution of the 
simulated sources, and (3) the presence of sources with 
absolute values of RMs larger than 300 rad m~^ which 
should be treated with suspicion due to the m: ambiguity 
in the definition of polarization angle. The total number 
of bins with errors increased by hand is about 10%. The 
algorithm is described in Appendix [Xj 

Some bins were completely removed from the analysis. 
First, these are the bins in polar regions with b > 50° 
(50 bins) and b < —40° (77 bins). In these regions the 
model values of RM are small because the projection 
of the MF onto the line of sight is suppressed (recall 
that in our model both disk and halo fields are paral- 
lel to the Galactic plane regardless of the parameters). 
On the other hand, there is strong evidence that some 
sort of 'screen' is pres ent in the direction of high lati- 
tudes (|Stil et al.l [20TTI ) which may provide a dominant 
contribution in these regions. The averaged RMs in po- 
lar regions, together with the assigned errors, are shown 
in Figures I13[ and [TJ] where some bias toward positive 
values of RMs can indeed be seen. 

Second, all bins with |6| < 10° were excluded because 
we could not trust the accuracy of the NVSS values in 
this region as explained in Section[2j Finally, we excluded 
all bins where the number of observed sources was smaller 
than 30 (this primarily concerns bins overlapping with 
the blind spot). Figure [5] shows the excluded bins in 
gray. . , , , 

A similar prescription was adopted for the analysis of 
the 1° strip along the Galactic plane using the RM set by 
KNMll, with the following changes. As the area in this 
case is 10 times smaller (~ 10 deg^ instead of 100 deg^) 
we set our lower limit to three sources in the bin. Also, we 
increased the assigned error value to ai instead of ai/3, 
thus attempting to reproduce the enhanced uncertainty 
in the disk due to the presence of a large number of local 
structures and very long propagation path length inside 
the disk (of order of tens of kpc). In the 1° strip, after 
visual examination, we removed 20 sources out of total 
number of 259. RM values of the removed sources clearly 
demonstrate large differences with RMs of neighboring 
sources. 



4.3. Fitting procedure 

The best values of the GMF parameters were identified 
by the maximum likelihood method. This amounts to 
minimizing the value of the defined as follows, 

2 _ (RMobs — R-Msim)? 

where A^tot is the total number of bins that we have in- 
cluded in the fit, (RMobs — RMsim)i is the difference be- 
tween the observed and simulated average RM values for 
the i— th bin and Si is the adopted error of the observed 
average. 

Since the RMs (Equation ([2|)) are hnear in the MF, 
the total RM of a given source is a simple sum of the 
disk and halo contributions. Making use of this prop- 
erty, we separately simulated more than 30,000 halo RM 
maps and 10,000 disk RM maps, corresponding to dif- 
ferent parameters varied within the ranges described in 
Section |3l At this stage the strengths of the disk and 
halo components were held fixed and set to values given 
in Section |3l Then each disk configuration was combined 
with each halo configuration with two arbitrary coeffi- 
cients a and /3. For a global GMF model thus obtained, 
the value of in Equation (fTTj) is a quadratic polyno- 
mial in a and /?, whose minimum with respect to a and 
/3 we found analytically. The resulting values of a and /3 
determined the strength of the disk and halo fields. 

The parameters of the GMF field were varied, within 
the ranges given in Section [31 by two different meth- 
ods: by uniformly covering a given parameter range 
with a regular grid (grid-type simulation) and by ran- 
domly picking parameter values from a uniform distri- 
bution over the desired range (Monte Carlo type sim- 
ulations). The disk field was simulated using the grid 
method, whereas for the halo we used both methods in 
roughly equal proportion. 

Because of the uncertainty in the absolute normaliza- 
tion of resulting from uncertainties in the error estima- 
tion (see Section H?2l) . we did not combine the Northern 
and Southern Galactic hemispheres in a single fit. The 
two hemispheres do not come on equal footing because 
of the blind spot in the Southern sky. Instead, our gen- 
eral fitting strategy was as follows. First we performed 
the fit using the NVSS data for the Northern hemisphere 
only, thus obtaining the preliminary set of models. Next, 
we fitted the NVSS data in the Southern hemisphere. 
We chose only the models that performed well in both 
tests and had the same value of the disk field. Note that 
this allowed for different halo fields in the Northern and 
Southern hemispheres. In these tests we imposed the 
additional constraint that |B| be greater than 1.8 /iG in 
both hemispheres, to conform with the results of other 
observations. Models were selected as acceptable if they 
satisfied the following requirement, 

xVXmin = Xrcd/X?od,min < 1 + V2/Nd.o.f., (12) 

where Xmin (x?od,min) ^ minimal value of (x?cd)- 
These are the models having parameters within ^ la 
from the best-fit ones. Note that the criterion is in- 
sensitive to the absolute normalization of the errors. We 
choose such a criterion to reduce the effect of the ambi- 
guities in the error estimation discussed in Section IT2l 
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Fig. 8. — Average rotation measures in bins. Red circles (blue 
squares) represent positive (negative) RMs. The color intensity 
reflects the absolute magnitude. Top: NVSS data. Bottom: best- 
fit model. 

Independently, we fitted the sources from the compila- 
tion by KNMll in the narrow strip of the galactic plane 
|6| < 1°. This imposed additional constraints on the 
structure of the field in the thin disk to which our fit of 
NVSS data is insensitive. Requiring that fits to both the 
NVSS and KNMll data be good, we performed the final 
selection of the models. 

Overall, we required that a successful model passes 
three separate fits: to the NVSS data in the Northern 
and Southern hemispheres, and to the KNMll data in 
the disk. Note that this is, in general, more restrictive 
than performing a single combined fit. 

5. RESULTS 

Figure[8]shows a comparison between one of the best fit 
models and the binned NVSS data. Positive average RMs 
are shown by red circles, and negative by blue squares. 
The intensity of the color indicates the absolute value of 
RM. It is seen that the general structure of the field is 
reproduced quite well. 

The results of the fits are as follows: 
North. The absolute minimum value of was equal to 
232 for 104 degrees of freedom (dof), x?od = 2.23. A 
bin-by-bin comparison between the data and the best-fit 
model is shown in Figure [S] Parameter ranges that were 
obtained from the fits are presented in Table 21 The fit 
was most sensitive to the pitch angle p. The vertical 
scale of the disk MF zq and the distance to the nearest 
reversal d affect the fit less. There was no significant 
infiuence of the Rc parameter, so we excluded it from 
further studies fixing its value at i?c = 5 kpc. The same 
concerns the parameters z^^^ and •^^2) "^^ich were fixed 



at z 



H 

1(1) 



— 0.25 kpc and z^^) = 0.4 kpc 



We could not make a successful fit using the disk com- 
ponent only. The addition of the halo component de- 
creases Xmin by a factor of 1.5. Due to a weaker effect 
of the halo MF, its parameters are estimated with larger 
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Fig. 9.— Best-fit model (blue solid line) vs. NVSS data (red 
points) in the Northern hemisphere. Bins not included in the fit 
are shaded with light blue. 

uncertainties than the disk ones. In addition, there is a 
degeneracy between the height of the halo z^ ^ the ampli- 
tude of the halo field and, to some extent, the radial size 
Rq . Even unrealistically large values of the halo field of 
order tens of fiG could give an acceptable fit if the halo is 
placed at a large distance from the galactic plane where 
the concentration of the electrons greatly decreases (cf. 
Figure [12]). There is also a degeneracy between the ver- 
tical scales of the GMF model (both halo and disk) and 
the vertical scale of the electron density (recall that the 
latter was held fixed during the fits). 

The strength of the disk field corresponding to the al- 
lowed set of models resides in the 1.8 2.8/iG range, 

largely concentrating around 2.0 /iG. The amplitude of 

the allowed halo field is in the range 2 12 ^G. Note 

that the halo field strength is constrained by the syn- 
chrotron radio-emission fro m the populati on of Galactic 
relativistic electrons (Hasl am et al.lTl982f ). These con- 
straints favor halo fields with the magnitude in the range 
2- -5 ^lG. 

Ring models of the disk field performed rather poorly 
in the Northern hemisphere no matter what was the 
halo field, giving the best value of Xmin = 288 for 101 
dof, which corresponds to Xred = 2.85 as compared to 
Xred = 2.23 for the best-fit model with the spiral disk. 
Thus, these models are disfavored by the NVSS data, 
confirming the earlier finding of Simard-Normandin & 
Kronbcrg(1980). 

Finally, our attempts to improve the quality of the fit 
by introducing a non-zero halo pitch angle were unsuc- 
cessful: there was no significant change of x^ when non- 
zero pitch angle of the halo field was varied within the 
range [—10°, 10°]. We also could not improve the fit by 
choosing the alternative model of the halo described by 
Equations ^ and (flOll . 

South. The data from the Southern hemisphere are 
more sparse because of the blind spot, and show more 
scatter, so the fits are generally worse. The best-fit ASS- 
type model gives x^ = 243 for 65 dof, which corresponds 



to x: 



red 



3.74. The resuhs for the best-fit BSS model 
are x^ = 235 for 65 dof (x?ed = 3.62). The bin-by-bin 
comparison between the data and the best-fit model in 
the Southern hemisphere is shown in Figure [TUl 

Disk. Finally, we fitted the RMs from KNMll at very 
low latitudes |6| < 1°. The best-fit ASS-type model gives 
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Fig. 10. — Best-fit model (blue solid line) vs. NVSS data (red 
points) in the Southern hemisphere. Bins not included in the fit 
are shaded with light blue. 
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Fig. 11. — Comparison between the best- fit model and the data 
from KNMll in the strip —1° < 6 < 1°. Areas not included into 
the fit are shaded with light-blue. 

= 30 for 29 dof (x?ed = 1-03), while the best-fit BSS 
model gives ^ = 37 for 29 dof (x?ed = 1-28). The 
comparison between these data and the best-fit model 
is shown in Figure [TT] The results are fully compati- 
ble with the other fits. They also a gree w ith the re- 
sults of iKronberg fc Newton-McGed (|2011[ ). Unfortu- 
nately, there is not enough information about the region 
around the direction to the Galactic center where the dis- 
crepancy between predictions of different models is the 
largest. This reduces the discriminative power of the 
fit. Nevertheless, the fit was very sensitive to the field 
strength in the disk, so its results impose most strin- 
gent constraints on this parameter. Also, it was shown 



by IKronberg fc Newton-McGed (|2011[ ) that the KNMll 
data favor the pitch angle near the disk of — 5.°5. Taking 
this constraint into account allows us to exclude positive 
pitch angles, which is not possible from the fits to the 
NVSS data alone. 

In both north and south fits of the NVSS data, the 
reduced values significantly exceed unity (especially in 
the Southern hemisphere). This could result either from 
our underestimation of errors, or from features in the 
data that cannot be accounted for by our models. Visual 
examination of the fits (Figures IH] and [TU| shows that 
there exists a systematic discrepancy between the data 

and the best-fit models in the region I « 130° 220° 

which corresponds to directions away from the Galactic 
center. The discrepancy is present in both Northern and 
Southern hemispheres. 

In fact, the observed variation of RM in the outer 
Galaxy is too strong to be accounted for in any model 
that has a thick disk, or halo with a nearly tangential field 
in the outer Galaxy (i.e., small pitch angle) if the field 
magnitude decays away from the Galactic center. So, ei- 
ther these models are completely wrong in this region, 
or the large RMs are caused by some nearby anomaly 
in the MF or the electron density, which is sufficiently 
close to subtend a large angle of about ~ 90°. The lat- 
ter option is supported by the fact that the variation of 
RM which is responsible for this feature does not visibly 
decrease away from the Galactic plane, as it would if it 
were due to an additional galactic arm or another global 
(and more distant) feature. Interestingly, the anomaly 
in the nearby re gion (the so-called region A) was point ed 
out previously ()Simar d-Normandin fc Kronbereill980[ ). 

For an additional check of this assumption we made a 
separate fit of RMs of the sources located in the direction 
of the inner Galaxy, i.e., at —90° < / < 90°. Excluding 
this apparent anomaly should improve the fit and show 
the extent to which the best-fit parameters were affected 
by its presence. The results of this additional test are as 
follows. 

North-inner. The best fit gives Xmin = 43 for 46 dof, 
considerably better than for the full sky. The allowed val- 
ues of parameters obtained from the restricted fit largely 
coincide with those obtained from the full fit. This indi- 
cates that the full fit is not dominated by the anomalous 
region. 

South-inner. Xmin = 35 for 22 dof. The quality of the 
fit has considerably increased as compared to the fit in 
the full range of longitudes, but is still worse than in the 
Northern hemisphere. 

All our attempts to fit the outer Galaxy separately 
(90° < / < 270°), even those including an extra disk 
component, were unsuccessful in the sense that we could 
not get Xrcd smaller than five. This confirms our conclu- 
sion that it is not possible to fit the anomalous region 
with the GMF models of the type considered here. 

The results of individual fits are presented in Table 2] 
in the form of acceptable ranges of model parameters. 
It should be stressed that there are dependencies and 
degeneracies among them; Table |4] only shows the max- 
imum variation of each parameter within the whole ac- 
ceptable set of models. In other words, not every combi- 
nation of the acceptable parameters gives an acceptable 
model. Parameters of the models which fit the data best 
and are in agreement with all fits are summarized in Ta- 
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TABLE 1 





P 
dcg 


20 

kpc 


d 
kpc 






ASS [-5, 


-4] U [3, 7] 


[0.7, 1.2] 


[-0.7,- 


■0.4] 


[1.8,2.2] 


BSS [-7, 


-3] U [5,8] 


[0.8, 1.2] 


[-0.9,- 


■0.5] 


[1.8,2.4] 



Note. — Ranges of parameters of the disk field corre- 
sponding to ~ la deviation from the best fit to the NVSS 
data and compatible with the KNMll data in the Galac- 
tic disk. Note that positive ranges of the pitch angle, al- 
though compatible with our fits, are strongly disfavored by 
the "shift-and-reflection" arguments of KNMll. 



TABLE 2 







kpc 


kpc 


RiT 
flG 








North 




vs. 


ASS 


[1.2, 2.4] 


[6.0, 15] 


[4.0, 12] 


vs. 


BSS 


[1.0, 2.0] 


[8.0, 15] 


[3.0, 12] 








South 




vs. 


ASS 


[1.0, 2.5] 


[3.5, 15] 


[1.5, 6.0] 


vs. 


BSS 


[1.0, 2.2] 


[4.0, 15] 


[2.5, 5.0] 



Note. — The same as Table [T] but for 
the halo parameters. 

blesIH and [21 

Finally, let us comment briefly on the other tests that 
we have made. First, as has been already mentioned in 
Section 13.31 we have run preliminary simulations with 
the original, lower vertical scale of electron density. For 
the realistic strength of the disk field Bq > 1.8 /j.G and 
the original electron density, we found only values 
larger than 390 as compared to 232 for the best-fit case. 
If we vary the magnitude of the disk field, the value of 
becomes close to the best-fit value for field strengths 
~ 0.7 /xG, which is clearly unacceptable. In conclusion, 
our simulations disfavor the original vertical scale of the 
electron density. 

Second, we have tested different combinations of sym- 
metries (with respect to the Galactic plane) of the halo 
and disk fields in order to check whether the NVSS data 
alone favor the symmetric disk and antisymmetric halo 
as we have chosen in our model. With our fitting pro- 
cedure which treats the Northern and Southern hemi- 
spheres independently, no additional simulations are re- 
quired for this test. For instance, for an anti-symmetric 
disk field we simply have to require that the coefficients 
with which the disk contribution enters the Northern and 
Southern hemispheres are opposite in sign. We could 
not find a satisfactory fit neither with ASS nor with BSS 
anti-symmetric disk field, obtaining the values of more 
than four times larger than the best-fit value. In the case 
of the symmetric disk and symmetric halo, the best fit 
gives 2.5 times larger than the best-fit model. Thus, 
the combination of the symmetric (with respect to b) disk 
and antisymmetric (with respect to b) halo is favored by 
the NVSS data. 

6. CONCLUSIONS 



The latest Faraday rotation data show an unmistakable 
coherent component of the RJMs over the sky. The KMs 
are correlated in sign and value over the largest angular 
scales. JVIoreover, the distribution of RlVIs is correlated 
with the Galactic plane and the direction to the Galactic 
center. This strongly indicates the Galactic origin of this 
coherence and provides firm evidence for the existence of 
a regular GIVEF. 

Apart from the Galactic plane region |6| < 10°, the 
global pattern of the RlVIs has different symmetry with 
respect to the Galactic plane in the inner and outer 
Galaxy: symmetric outside and anti-symmetric inside 
(see Figure [5]). This change of the symmetry cannot 
be achieved with a single (symmetric or antisymmetric) 
component of the GIVIF. At least two components — disk 
and halo — are required (Figure |4]). Since the study of 
the MF in the disk shows that it does not change sign 
across the Galactic plane itself, the combination of the 
symmetric disk field and the anti-symmetric halo field 
is favored. We have demonstrated that such a combi- 
nation, with properly chosen amplitudes of components, 
reproduces the observed pattern of RJVIs very well. 

Several further conclusions can be made. 

1 . The cross check between the NVSS and the smaller 
KNJVIll set shows that the quahty of the NVSS is 
sufficient for our purposes all over the sky except, 
in the |5| < 10° band around the Galactic plane. 

2. There is a feature in the observed RlVIs near the 
direction to the Galactic anti-center that cannot 
be accounted for by models of the type considered 
in this paper (nearly tangential field decaying away 
from the Galactic center). This feature may be due 
to some unspecified local structure. 

3. A ring model of the GMF is disfavored as compared 
to the model with logarithmic spiral arms. In the 
former case the best value of Xrcd ^ factor of 
> 1.3 larger than in the latter. 

4. Our fits favor the increased vertical scale of 
the electron d e nsity distribution as proposed by 
IGaensler et aI1(|20?)8h . We were unable to find any 
fit with the original vertical scale of NE2001 model 
and stron g disk field. This is i n agreement with the 
resuhs of lSun fc ReichI (IMol) . 

5. Our fits favor somewhat smaller values of the pitch 
angle p than have been discussed in the literature 
previously. Th ese smaller values agree with t he re - 
cent results of iKronberg fc Newton- IVTcGed (|2011l ) 
where the reference value of p = —5.° 5 was ob- 
tained by a totally different method. 

Our fitting procedure could not discriminate between 
the ASS and BSS models because a major contribution 
to the extragalactic source RJVIs (especially at \b\ > 15°) 
comes from the near vicinity of the solar system; in this 
region both models behave quite similarly. A more de- 
tailed analysis of the Galactic plane region may be more 
sensitive to the differences between the ASS and BSS 
models. 

Two "benchmark" models (one with the ASS and one 
with the BSS disk field) are presented in Table [3] These 
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TABLE 3 
Benchmark model parameters. 



ASS 



BSS 



Disk: 
P 

d 

Bo 

Rc 

Halo (North): 

% 
^i(i) 

1(2) 

Halo (South): 

1(2) 



-5° 
1.0 kpc 
-0.6 kpc 
2.0 

5.0 kpc 



1.3 kpc 

8 kpc 

0.25 kpc 
0.4 kpc 



1.3 kpc 
8 kpc 
2^jG 
0.25 kpc 
0.4 kpc 



-6° 
1.0 kpc 
-0.6 kpc 
2.0 AiG 
5.0 kpc 



1.3 kpc 
8 kpc 
4AiG 
0.25 kpc 

0.4 kpc 



1.3 kpc 
8 kpc 
4mG 
0.25 kpc 

0.4 kpc 



models fit the RM data well and have values of param- 
eters favored by other studies (not too small pitch angle 
p, not too strong halo field and the magnitude in the 
vicinity of the Earth close to 2 \iG). They can be used, 
e.g., for propagation of cosmic rays through the Galaxy. 
At ultra-high energies, an accurate three-dimensional de- 
scription of the Galactic magnetic field may be crucial for 
identification of extragalactic cosmic ray sources. 

In this paper we used only the Faraday RMs of ex- 
tragalactic radio sources. More precise understanding of 
the GMF can be achieved by taking into account of other 
data, also sensitive to the MF, such as RMs of pulsars, 
measurements of the galactic synchrotron radiation, etc. 
This would permit better constraints on the unknown 
quantities such as the three-dimensional electron density 
profiles and other properties of the ISM. It will also help 
to resolve some degeneracies inherent in our models. 
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APPENDIX 



ERROR ASSIGNMENT. EXCLUDED BINS AND BINS WITH MANUALLY ASSIGNED ERRORS 

In our fitting procedure, each bin included in the analysis was characterized by an average RM and its error (5. Here 
we describe the algorithm by which these errors were assigned on a bin-by-bin basis. 

First, we estimated the intrinsic dispersion of RMs of extragalactic sources. To this end we investigated the regions 
on the sky where the influence of the GMF is expected to be suppressed. Since RMs are proportional to the MF 
component that is parallel to the line of sight, these are the regions around the Galactic poles. Considering the sources 
with |6| > 80° we found that the dispersion of RMs of individual sources in these regions was about CTqso = 15 radm~^. 
Comparing to Figure [3] we see that the dominant contribution to this number comes from the measurement errors. In 
the error assignment algorithm, we adopted the value of CTqso as an indicator of whether a given bin is perturbed by 
effects other than the random scatter of sources. 

If a bin showed the dispersion of RMs Oi smaller than cTqso we considered it unperturbed and assigned its error on 
a statistical basis. Since the average of a sample of N elements varies around the true average with the width a/y/N ^ 
where is the variance of the distribution, we assigned to such bins an error of 6i = crqso/\/Ni, Ni being the number 
of sources in the bin. In practice, this situation occurs only for high-latitude bins. 

In the opposite case, i.e., when the dispersion ai was larger than Cqso, we considered the bin as perturbed by coherence 
effects. These may be due to, e.g., random MF with the coherence length comparable to the bin size, fluctuations of 
electron density, strong variation of the regular component across the bin, etc. In this case the statistical approach 
is not correct, and a precise determination of the error is difficult. We have assumed that the error in such bins 
is some fraction of the bin variance ai, which we have adopted to be 1/3, so that the error in perturbed bins was 
taken Si = ai/3. We cannot rigorously derive this coefficient (although it can be motivated by a toy model presented 
below). Since this uncertainty propagates directly into the value of x^i it results in some uncertainty of the absolute 
normalization of the latter. Note, however, that the ratios of 3,re less affected by this uncertainty. 

The following toy model shows how the proportionality between the error and the dispersion in the bin may come 
about. Let the random part of the MF consist of cells of a typical size d ^ 100 pc, such that their contributions to 
RM have zero mean and the dispersion ctq. If the cells were statistically independent, individual sources would have 
an RM dispersion given by ti = o'oVA'tot, where iVtot is the total number of cells along the line of sight (we have 
neglected the internal scatter). However, for sources in one bin a number A'coh of nearby cells contributes coherently. 
If A'coh ^ A'tot this will not change the dispersion a over the bin. The error, on the contrary, will be dominated by 
the coherent contribution, 6 ~ croV^coh = cr-\/7Vcoh/A^tot- Thus, the error and the dispersion are proportional to one 
another with the coefficient y/ Ncoh/Ntot- The latter coefficient is purely geometrical. Taking A'coh ~ 1/^*, where 9 
is the bin size, and A'tot ~ L/d, where L is the length of the trajectory inside GMF, we find S ~ cr/3. Interestingly, 
similar conclusions (proportionality and a similar coefficient) are obtained by considering a perturbation caused by a 
close group of a few sources having similar RMs substantially different from the rest of the sources in the bin. 
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As a last step, all bins were subject to visual scrutiny, with the result that some of the bins were completely removed 
from the analysis. These are the bins showing obvious pathologies like small number of sources or non-uniform 
distribution across the bin (e.g., some bins at the boundary of the NVSS blind spot), compact aggregates of sources 
with RMs drastically different from the surrounding value, etc. Finally, some bins with such pathologies were assigned 
large errors (2-3 times larger than prescribed by the above algorithm) by hand. This was done in cases where the 
corresponding bin might serve as a good discriminator between models, even with a large error. A constant error of 
40radm~^ was used in this case. This value was chosen in such a way that these bins do not contribute significantly 
into the value of x^, at least near the best-fit point. We have checked that the results of the fits do not depend on the 
precise choice of these errors. The number of such bins is about 10%; they are easily identifiable in Figures [6] and [T] 

Here is the complete list of excluded bins and bins with manually assigned errors. 

Northern hemisphere. An error of 40radm~'^ was assigned to the bins in the region 30° < I < 60°, 10° < b < 40° 
(local structures described in (Wolle ben et al.ll2010D '): to the bins in the area -115° < / < -90°, 20° < < 30° (a 
prominent structure with large positive RMs, possibly related with the Gum nebula); to the bin located at (15°, 25°) 
(strong non-uniformity: only 12 sources are located at 10° < b < 15°, while 38 at 15° < b < 20°). Excluded from the 
analysis were the bins in the area —10° < I < 30°, 20° < b < 30° (unresolved structure in RM and a large positive 
spot at (-10°, 25°)); bins in the area 70° < I < 90°, 10° < 6 < 20° (large negative spot centered at (80°, 15°); bin 
affected by Gum nebula at {I = — 100°,& < 20°); the bin centered at (—30°, 15°) (strongly non-uniform distribution). 
Southern hemisphere. An error of 40 rad m~^ was assigned to the bin (105°, —15°) (local structure resulting in positive 
excess ); to the bin (165°, — 15°) (positive spot); to the bin centered at (180°, -15°) (negative spot near the anti-center 
direction); and to the bin (170°, —25°) (large positive spot). 
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TABLE 4 

Ranges of disk and halo parameters of the acceptable models 









P 




d 


OL 








/3 






deg 


kpc 


(kpc) 




kpc 


kpc 














North 










A 


240 


2.30 


[-10, -3] 


[0.5, 1.3] 


[-1.1, -0.3] 


[0.9, 2.0] 


[1.0, 1.4] 


[4, 15] 


[0.4, 0.7] 


B 


232 


2.23 


-fi -101 


[0.5, 1.3] 


[-1.1, -0.5] 


[0.9, 1.5] 


10 14 


6 1 5 


[0.3, 1.0] 












South 










A 


243 


3.74 


[-10, 10] 


[0.5, 1.5] 


[-1.1, 1.1] 


[0.9, 2.0] 


[1.0, 3.5] 


[3.5, 15] 


[0.2, 2.0] 


B 


235 


3.62 


[-10, 10] 


[0.5, 1.1] 


[-1.0, 1.0] 


2 2 


[1.0, 3.5] 


3 5 15 


2 2 












Plane 










A 


30 


1.03 


[-7, 10] 


[0.5, 1.4] 


[-1.4, 1.4] 


[0.9, 1.1] 








B 


37 


1.27 


Ls -'^lLJf5 Si 


[0.5, 1.5] 


[-0.8, 1.2] 


[1.0, 1.2] 


















North-inner 










A 


43 


0.93 


[-5, 6] 


[0.5, 1.2] 


[-1.0, 0.9] 


[0.9, 3.0] 


[1.2, 2.4] 


[6.0, 8.0] 


[0.4, 1.2] 


B 


45 


0.98 


[-7, 9] 


[0.5, 1.4] 


[-0.6, 1.1] 


[0.9, 1.1] 


[1.0, 2.0] 


[8.0, 15] 


[0.4-1.2] 












South-inner 










A 


40 


1.81 


[-7, 11] 


[0.5, 1.3] 


[-0.2, -0.9] 


[0.9, 3.5] 


[1.0, 2.5] 


[8.5, 15] 


[0.1, 0.6] 


B 


35 


1.59 


[-5, 9] 


[0.5, 1.2] 


[-0.3, -0.8] 


[0.9, 1.2] 


[1.2, 2.2] 


[3.5, 5.5] 


[0.2, 0.5] 



Note. — A and B denote ASS and BSS models, correspondingly. The parameters a and /3 determine 
the strengths of the disk and halo components as explained in Sect. 14.31 




5 10 15 20 

B, IIG 



Fig. 12. — Height of the halo Zq vs. its strength Bq . Each dot represents an acceptable model. The degeneracy could be easily seen: a 
very strong halo field gives an acceptable fit provided the vertical scale is high. 
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Fig. 13. — NVSS data for b > 50°. The bias towards positive values of RMs, perhaps due to the small vertical component Bz, could be 
seen. 
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Fig. 14. — NVSS data for b < —50°. Bias towards positive values of RMs could be seen here as well. 

REFERENCES 



Beck, R. 2001, Space Sci. Rev., 99, 243 

Beck, R. 2008, in American Institute of Physics Conference 
Series, Vol. 1085, American Institute of Physics Conference 
Series, ed. F. A. Aharonian, W. Hofmann, & F. Rieger, 83-96 

Brandenburg, A., Donner, K. J., Moss, D., et al. 1992, A&A, 259, 
453 

Braun, R., Heald, C, & Beck, R. 2010, A&A, 514, A42+ 
Brown, J. C, Haverkorn, M., Gaensler, B. M., et al. 2007, ApJ, 
663, 258 

Brown, J. C, Taylor, A. R., & Jackel, B. J. 2003, ApJS, 145, 213 
Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 
115, 1693 

Cordes, J. M., &; Lazio, T. J. W. 2002, |arXiv:astro-ph/02 07156| 
Crutcher, R. M. 1999, ApJ, 520, 706 

Frick, P., Stepanov, R., Shukurov, A., & Sokoloff, D. 2001, 

MNRAS, 325, 649 
Gaensler, B. M., Madsen, G. J., Chatterjee, S., &; Mao, S. A. 

2008, Publ. Astr. Soc. Australia, 25, 184 
Han, J. L., Manchester, R. N., Berkhuijsen, E. M., & Beck, R. 

1997, A&A, 322, 98 
Han, J. L., Manchester, R. N., Lyne, A. G., Qiao, G. J., & van 

Straten, W. 2006, ApJ, 642, 868 
Han, J. L., & Qiao, G. J. 1994, A&A, 288, 759 
Haslam, C. G. T., Salter, C. J., Stoffel, H., & Wilson, W. E. 1982, 

A&AS. 47, 1 



Heiles, C. 1996, ApJ, 642, 316 

Jaffe, T. R., Leahy, J. P., Banday, A. J., et al. 2010, MNRAS, 
401, 1013 

Jansson, R., Farrar, G. R., Waelkens, A. H., & Enfilin, T. A. 

2009, JCAP, 7, 21 
Klein, U., Mack, K., Gregorini, L., & Vigotti, M. 2003, A&A, 

406, 579 

Krause, M. 2009, in Revista Mexicana de Astronomia y 

Astrofisica Conference Series, Vol. 36, Revista Mexicana de 
Astronomia y Astrofisica Conference Series, 25-29 

Kronberg, P. P., & Newton-McGee, K. J. 2011, Publ. 
Astr. Soc. Australia, 82, 171 

Mao, S. A., Gaensler, B. M., Haverkorn, M., et al. 2010, ApJ, 
714, 1170 

Mao, S. A., Gaensler, B. M., Stanimirovic, S., et al. 2008, ApJ, 
688, 1029 

Men, H., Ferriere, K., & Han, J. L. 2008, A&A, 486, 819 

Nishiyama, S., et al. 2010, ApJ, 722, L23 

Prouza, M., & Smi'da, R. 2003, A&A, 410, 1 

Ruiz-Granados, B., Rubino-Martm, J. A., & Battaner, E. 2010, 

A&A, 522, A73+ 
Simard-Normandin, M., & Kronberg, P. P. 1979, Nature, 279, 115 
— . 1980, ApJ, 242, 74 

Simard-Normandin, M., Kronberg, P. P., & Button, S. 1981, 
ApJS, 45, 97 



Deriving structure of the Galactic Magnetic Field 



15 



Sofuc, Y., & Fujimoto, M. 1983, ApJ, 265, 722 
Stanev, T. 1997, ApJ, 479, 290 

Stil, J. M., Taylor, A. R., & Sunstrum, C. 2011, ApJ, 726, 4 
Sun, X., & Reich, W. 2010, Research in Astronomy and 

Astrophysics, 10, 1287 
Sun, X. H., Reich, W., Waelkens, A., & EniSUn, T. A. 2008, A&A, 

477, 573 

Taylor, A. R., Stil, J. M., & Sunstrum, C. 2009, ApJ, 702, 1230 



Tinyakov, P. G., & Tkachcv, I. I. 2002, Astroparticle Physics, 18, 
165 

— . 2005, Astroparticle Physics, 24, 32 
Vallee, J. P. 2005, ApJ, 619, 297 

Wielebinski, R. 2005, in Lecture Notes in Physics, Berlin Springer 
Verlag, Vol. 664, Cosmic Magnetic Fields, ed. R. Wielebinski & 
R. Beck, 89-+ 

WoUeben, M., et al. 2010, ApJ, 724, L48 



